Sparsity-Based Data Reconstruction Models for Biomedical Imaging

We propose new regularization models to solve inverse problems encountered in biomedical imaging applications. In formulating mathematical schemes, we base our approach on the sparse signal processing principles that have emerged as a central paradigm in the field. We adopt a variational perspective and specify the proposed sparsity-promoting data reconstruction models as energy minimization problems. To design practical algorithms, we develop novel iterative methods to efficiently perform the consequent optimization task. The thesis is organized in three main parts. In the first part, our main contribution is the introduction of a proper statistically-based discretization paradigm for inverse problems. In particular, our framework considers a continuous-domain stochastic signal model and characterizes a specific class of inference algorithms. We show that derived inference-based methods cover the classical Tikhonov-type techniques as well as a wide range of the sparsity-promoting schemes including the well-known `1-norm regularization and its nonconvex variants. This provides a unifying stochastic perception of the resolution of inverse problems. In the second part, we propose a novel phase retrieval algorithm for imaging unstained biological samples that are optically thin. In specific terms, we use the transport-of-intensity equation (TIE) relating the spatial phase map of a field to the derivative of its intensity along the propagation direction. We analyze the implications of using the TIE formalism with finer and coarser approximations of said derivative. Based on this analysis, our contribution is a practical phase reconstruction algorithm that incorporates a sparsity-based regularization. The developed technique operates with a standard bright-field microscope. Experiments on real data illustrate that our phase reconstruction algorithm is viable and can be a low-cost alternative to dedicated phase microscopes.